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The sit(2|l) coherent-state path-integral representation of the partition function of the t — J 
model of strongly correlated electrons is derived at finite doping. The emergent effective action is 
compared to the one proposed earlier on phenomenological grounds by Shankar to describe holes 
in an antiferromagnet (Nucl.Phys. B330 (1990) 433). The t — J model effective action is found 
to have an important "extra" factor with no analogue in Shankar's action. It represents the local 
constraint of no double electron occupancy and reflects the rearrangement of the underlying phase- 
^ ■ space manifold due to the presence of strong electron correlation. This important ingredient is 

shown to be essential to describe the physics of strongly correlated electron systems. 
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■ I. INTRODUCTION 

In this work, we discuss a path-integral representation of the partition function for strongly correlated electron 
^ I systems. In particular, we are interested in the effective low-energy action to describe a lightly doped antiferromag- 
net (AF). Strong electron correlations essentially determine the low-energy physics of high-Tc superconductors |2|, 
heavy-fermion systems 0], itinerant ferromagnets [1], as well as some optical lattices @. Accordingly, many various 
approaches have been proposed to study those corresponding low-energy effective actions. We will comment on some 
, of them at the end of the paper. 

■ Our work has actually been motivated by Shankar's contribution 1]. Namely, two decades ago Shankar put forward 
I , a conjecture assuming that the low-energy action to describe a doped quantum AF involves spinless fermions locally 

' coupled to a compact U{1) gauge field. This gauge field is driven by the AF fluctuations controlled by a nonlinear 
] sigma-model. This approach was further discussed in Q . A natural question then arises as to whether that action can 
O . be derived directly from the t — J model for strongly correlated electrons at finite doping. After all, this microscopical 
model is widely believed to capture the essential physics of a lightly doped AF. 

In the present paper we show that Shankar's effective action can indeed be derived from the su(2|l) path-integral 
representation of the partition function of the t — J model. However, the resulting path-integral measures differ by 
an important factor. This distinction is a manifestation of the strong coupling nature of the problem due to the 
no double occupancy (NDO) constraint. It may seem that this constraint is already fully accounted in Shankar's 
I _ theory by the requirement that the fermions are spinless, since no double fermion occupation is possible in this case, 
ly-j ■ However, the NDO constraint modifies the original on-site Hilbert space as well as the on-site path-integral phase 
space, which is not explicitly taken into consideration in [Ij. Shankar proposed instead an effective action right in 
terms of conventional fermion and spin fields to describe a doped AF. This of course implies a standard measure in 
the path integral. Although, his action may presumably describe some unconstrained spin-fermion models, it is not 
appropriate for strongly correlated electrons. 

Our aim is to demonstrate in what way the NDO constraint modifies the theory discussed in [1]. To start with, 
we illustrate our point by considering a path-integral representation of the partition function for a simple single-site 
Hamiltonian. After that, we address the t — J model close to half filling which is precisely the physically most relevant 
' situation for strongly correlated electrons. This model is believed to capture the low-energy physics of lightly doped 
quantum AF. 

Quite plainly, an electron system is said to be strongly correlated if the leading energy scale in the problem is the 
on-site Coulomb repulsion energy U . In this case the low-energy sector of the underlying on-site Hilbert space should 
be modified to exclude doubly occupied states. Such a modification results in an entirely new physics to account for 
the relevant low-energy excitations. Formally, strong correlations are encoded into the projected electron (Hubbard) 
operators. They act directly in the restricted Hilbert space as opposed to the conventional electron operators which 
describe the unconstrained system. In contrast with the conventional fermion operators which generate the standard 
fermionic algebra, the new operators obey more complicated commutation/anticommutation relations and are closed 
into a superalgebra su{2\\). It is therefore natural to seek a path-integral formalism that takes the structure of that 
superalgebra fully into consideration. In analogy with the conventional su{2) spin path integral Q, this can be carried 
out by employing an appropriate coherent-state basis associated with the sm(2|1) superalgebra representations. To 
make our exposition self-contained we employ some notation and results in connection with a continuum sm(2|1) 



2 



path integral already reported elsewhere [8|. However, the conclusion we reach below requires a more sophisticated 
approach based on a carefully defined time-lattice representation of the corresponding path integral. 

Just to get an idea about what strong correlations really are, consider the U — oo Hubbard model which is known 
to capture the extreme limit for strongly correlated electrons. The Hamiltonian reads 

ij^a i cr 

Here c^j^^{ci„) is an on-site creation (annihilation) operator of an electron excitation with the spin projection a =t,i- 
The hopping amplitudes tij represent jumps between nearest-neighbor (nn) and next-to-nearest-neighbor (nnn)sites 
and are zero otherwise. Throughout this paper we will be considering a D dimensional bipartite lattice, L = A® B. 
The chemical potential /i controls the total number of vacancies (empty sites). 

The infinitely large on-site Coulomb repulsion is accounted for by the local NDO constraint, < 1, that restricts the 
on-site Hilbert space to states with at most one electron per site. This is the essence of strong electron correlations. 
In the absence of this constraint, this model simply describes a system of noninteracting electrons and it can be 
trivially diagonalized in the momentum space. The NDO constraint makes this problem nontrivial and its exact 
solution is still unknown for spatial dimensions > II?. The physics behind the model (HI) is certainly far from trivial. 
Indeed, one of the few exact results was proved by Nagaoka Q who showed that for one hole the ground state of the 
U = oo Hubbard model is a fully saturated ferromagnet. This provides an interesting example of a quantum system 
with ferromagnet ordering due to a purely kinetic-energy effect driven by hole hopping (itinerant ferromagnetism). 
Unfortunately, despite very extensive work over the years, both this model and itinerant ferromagnetism are still 
poorly understood. One of the open important questions related to this concerns the thermodynamic stability of 
the Nagaoka phase. That is, whether or not the Nagaoka state is stable when the density of holes is finite in the 
thermodynamic limit. 

The local NDO constraint can be explicitly incorporated into the theory by projecting the Hamiltonian ([l} onto 
the restricted Hilbert space J^^ , where the 3D on-site Hilbert space 7^^''^* is spanned by the vectors |0)i (empty 

site), I t)i (spin-up electron), and | (spin-down electron): 



H ^ H = VHV, r = Y[ri 



where the Gutzwiller projection operator Vi is given by 
so that 

In this 3D subspace the constrained electron operators can be identified with the Hubbard operators, 

xo- = (x-o)t = |o)(a|, <7=t,;. 

It follows that Cio- = X^"' , provided rii^rii^^ = 0. The last requirement eliminates the doubly occupied fermionic 
states. As a result the Hamiltonian ([T]) in the constrained physical space takes the form 



= ^tY^XfXf -i^Y.^^- (2) 

zj.cr i 

Here X"" = X^" X"^ — |0) (0| = 1 — X^o- ^a^<y stands for the on-site vacancy number operator, so that the concentration 
of vacancies becomes ^ = ^ 'Yliik^f^)- The on-site operator X^" removes an electron with the spin projection a and 
creates a vacancy (empty state) which is a spin singlet. 

The important point is that the fermionic Hubbard operators X'^'^ = {X"''^y along with the bosonic ones, X"'"' , X*"^, 
are closed under commutation/anticommutation relations into the superalgebra su(2|l) The su(2|l) superalgebra 
can be thought of as the simplest possible extension of the conventional spin su{2) algebra to incorporate fermionic 
degrees of freedom. Namely, the bosonic sector of the su(2|l) consists of three bosonic superspin operators, 

Q+ = X^\ Q-=X^\ g-^ l(Xtt_x^i) (3) 
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closed into su{2), and a bosonic operator that generates a factor of the maximal even subalgebra su{2) x u{l) 
of sm(2|1). The fermionic sector is constructed out of four operators X"''^, that transform in a spinor representation 
of su{2). We do not intend here to discuss a general theory of the su(2|l) irreducible representations that can be 
found elsewhere We will focus instead on a specific lowest 3_D irreducible representation of sm(2|1) generated 

by the Hubbard operators. This representation acts in the physical on-site Hilbert space of model ([ij with the basis 
vectors, | t)i I i): |0), and consists of nine operators X^^ — | A) (A'|, A, A' = 0, a. The resolution of unity in this space, 

/ = io)(oi + it>(ti + u)ai = E^''' (4) 

A 

singles out eight independent generators of su{2\l) listed above. The special important property of this representation 
is that X^^'X^"^'" = 6x'\"X^^"' . 



II. s«(2 1) COHERENT-STATE MANIFOLD 

The normalizable coherent states (CS's) associated with the lowest irreducible representation of sw(2|l) superalgebra 
spanned by Hubbard operators take the form 

\z,0 = (l + zz + eO"'/'exp(zX^t + ^XOt)|t) 

= (i + zz + ee)-'/'(it> + ^i;> + eio)), (5) 

where a complex even Grassmann parameter z and an odd complex Grassmann parameter ^ are the inhomogeneous 
(proper) coordinated of a point on a supersphere, {z,^) g S*^'^ ~ CP^'^ — SU{2\1)/U{1\1). Here CP^'^ stands for 
a complex projective superspace with a complex dimension (1, 1). It can be thought of as a minimal superextension 
of an ordinary projective space CP^ homeomorphic to a two-sphere, CP^ ~ S^. The odd Grassmann parameter 
appears in ([S]) due to the fact that X^^ is a fermionic operator in contrast with the operator X^"^. The product 
represents therefore a bosonic quantity as required. 

One can use as well the homogeneous coordinates on the supersphere, so that z = z^/z"^,^ = 0/z'^, with 

z'^ ^ 0. Then the coset CP^'^ manifold is defined by the equation 

\z^\^ + \z^\^ + 99 ^1. 

The supergroup SU{2\1) acts on CP^I^ according to Z ^ ^ gZ, whereZ = {z^,z'^,9y and g G SU{2\1). This 
generates a corresponding transformation of the inhomogeneous coordinates, {z,£,). In particular, if one chooses g to 
represent a pure spin rotation, 

7i\y ("-)^™'^'- 

one gets 

uz + v i 
—vz + u —vz -t- u 

Note that both the bosonic and fermionic fields transform themselves under SU (2) spin rotations. 
At ^ = 0, the su(2|l) CS reduces to the ordinary spin su(2) CS, 

\z, e = 0) = ^ \z) = -7=^=rm ^Mzs-)\ t> - J-Ti2 ^\ ^> + ^1 

where the complex number z is a stereographic coordinate of a point on an ordinary sphere, z G 5^ ~ CP^ = 
SU{2)/U{1). The spin operators S obey the standard commutation relations 

[S,,S±] = ±S±, [^+,5_] = 25„ 5^2^3/4. (8) 

These operators coincide with bosonic generators Q, of su(2|l) at half filling, in which case the on-site Hilbert space 
is reduced and spanned only by the vectors | t), | |). In contrast, at z = 0, the state |^) = |z = 0, ^) represents a pure 
fermionic CS. 
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At this stage it is helpful to introduce the important notion of the covariant (Berezin) symbol for a Hubbard 
operator X. It can also be referred to as a coherent-state symbol and is defined as follows 

X,o.:^{zX\X\z,Ci. (9) 

Explicitly, we find 

l+|z|2' l+|z|2' 



coi; 1(119' cov 1 I I 12 ' \ ^/ 



where the covariant symbol of the hole number operator reads 



The corresponding CS symbols of the su(2) generators are evaluated to be {Scov ■— (^l-S'lz)): 

S+ - ' S- - ^ =lf}—K.\ (13) 

The important point is that the covariant symbol of the Hamiltonian associated with the algebra generators enters 
the path-integral action for the partition function. For a compact simple algebra those symbols are in one-to-one 
correspondence with the algebra generators. This is the case for both the sm(2|1) and su{2) algebras. Note also that 
at half filling, S = 0, we get = so that the symbol of the Q operator reduces to that of the conventional spin 
operator, Qcov = Scov In general this is not the case. Although the generalized spin operators Q fulfill the same 
commutation relations as operators S do, the operator is not a c-number. Namely, = S^{1 -t- X^"). This means 
that operators Q are closed into an algebra larger than su{2). 

The su(2|l) CS is parameterized by the coordinates of a point on the coherent-state supermanifold CP^'^. The 
latter appears as a classical phase space for the constrained electrons. The symplectic structure is given by a canonical 
symplectic two-form, 

n = d{i{z, ^\d\z, 0) = -i [dz{a)dz + dz{a)d^ + d£,{j)dz + d^{y)d^ , (14) 
where the external derivative 

J ^9 ^ d -d 

oz oz d£, Ot, 

and 



(l+|^|2+^-^)2' (l+|^|2+^-^)2' ' (l + |^|2+^-^)2'« 1+|^|2- 

We use the convention that a superform on CP^'^ is Z x Z2 graded, where the Z- gradation is the usual gradation of 
the de Rham complexes, while the Z2-gradation is a natural gradation of Grassmann algebra [ll|. We thus have 

dzdz = —dzdz, dzd^ = —d^^dz, d^^d^ — d^d(^. 

Classical dynamics on CP^'^ is governed by the Poisson brackets generated by the symplectic structure. Namely, 
for any two functions g and h that represent classical observables on CP^'^, we get 

{g,h}pB = +iXg\dh, (15) 

where Xg denotes a vector field associated to g and the symbol J stands for the interior product (or contraction), so 
that [12] 



Xg\n = dg. 
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In particular, the covariant symbols (|10lll|) are closed into the sm(2|1) superalgebra under the action of brackets (|15l) . 
Up to a multiplicative constant the SU{2\1) invariant measure is given by 

.. .. - dzdzd^dS, , , 

d^i cx sdet\\il^i\\dzdzd^d^ cx — — — — (16) 

1 + zz + 44 

where "sdet^^ stands for the superdeterminant (or Berezian) jT3| . while stands for the supermatrix form of £7, 

namely — dx^^l^gda;^, x°- — z,£_. 



III. su{2\l) PATH INTEGRAL: PRELIMINARIES 

We now seek a path-integral representation of the partition function Z = tre~^^ where H is expressed in terms of 
Hubbard operators. This implies that the NDO constraint is explicitly resolved from the outset, so that no redundant 
gauge dependent variables emerge within that approach. To work out the path-integral formalism to deal directly with 
the Hubbard operators, consider the on-site 31? physical Hilbert space Hphys spanned by the vectors |A), A —'f,l,0. 
The Hubbard operators act in this space and are naturally Z2 graded. As was already mentioned, there are bosonic 
or even-graded and fermionic or odd-graded operators. Accordingly, the basis in Uphys becomes graded as well. We 
choose the spin-up | f) and spin-down | I) states to be even-graded, whereas the empty state (vacancy) |0) is considered 
to be odd-graded. This means that odd Grassmann parameters commute with the spin states and anti-commute with 
the vacancy vector. From now on we consider the grading of the physical 3Z) subspace to be fixed in accordance with 
this. For any operator H acting in Hphys one then gets 

tre-^" = J2{X\e-^"\X), (17) 

A 

and 

stre-^" = ^(-)9™'*(l^>)(A|e-^^|A), (18) 

A 

By using the su(2|l) CS (O, Eg. p?)) can be rewritten in the form 
Here the S'?7(2|l) invariant measure 

_ dzdzd^d^ 
dtisu{2\i) =d^i^ - (20) 

27ri(l + zz + 44) 

agrees with Eq. ()16|) . Integration over z,z is performed on a complex plane. The integration over Grassmann 
parameters is carried out according to the standard rules [l3|. 



^ - / = 0, / dee = / = 1. 

Notice that in contrast with Eq. ([T9|) 

stre-^" ^ I d^l{z,C\e-^"\z,0. (21) 



The distinction between Eqs. (|19|) and pip results in different boundary conditions (anti-periodic versus periodic) 
for Grassmann amplitudes in the corresponding path integrals for the partition/superpartition functions. 

In the course of the derivation of the partition function we essentially rely on a certain change of the path-integral 
variables. To justify it, we start with an example of a simple single-site Hamiltonian that also admits a treatment 
in terms of ordinary multiple integrals. After that, we reconsider the problem in terms of the time-continuous path 
integral. 

Consider the single-site Hamiltonian 

H = tiX°°. (22) 
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Operator = |0)(0| is represented in Hphys by 3 x 3 matrix with eigenvalues 0, 0, 1. As a result, the partition 
function reduces to 



and the vacancy occupation number becomes 



Zo^2 + e 



1 



1 



This is to be compared with the conventional spinless fermion occupation number, 

1 



(23) 



(24) 



e^A' + 1 ■ 



Since (X"")"^ = X™ we get 
Equation (fTH)) then gives 

Here we have used that 



= l + X°"(e-^'^-l 



dM(2,e|(l + X00(e-^^-l))|z,-e) 



dzdzd^d^ 



27ri(l + ZZ + CO 



(e-*^^ - 1) 



1 + zz + CC 1 + zz + CC 



l + z'z + 



(25) 



(26) 



^/1 + z'z' +ec'V'^ + zz + cc 

Although integral ([25]) can be calculated immediately, it is instructive at this stage to make a change of variables 

z^z, C^^v/l + I^P (27) 
to bring the measure into a more tractable form. 



dzdzd^d^ 



dzdzdS^dS, 

27ri(l + zz + f^) ' 2iTi{l + zz)' 



(28) 



Equation ((25)) then becomes 



dzdzd^d^ 
27rj(l + ^z)2 

dzdzd^d^ 
27ri(l + zz)2 



(1-^0 



ee(e-^^ - 1) 



[l-3^^-^"e(e-^''-l)]) = 3 + . 



- 1 = 2 + e" 



(29) 



which agrees with p3|) . 

Under the change (|27|) we get X^^^ — >• so that it may seem we have reduced the problem to that of spinless 
fermions. However, this is not the case. The measure (|28p involves the extra factor (1 — ^^) comparing with the 
standard product of the SU{2) invariant spin and spinless fermion measures. 



dzdz 



27ri(l + zz)' 



;d^dC 



This extra factor reflects a nontrivial symplectic structure of the underlying path-integral phase space. Namely, the 
transformation (j27p brings the symplectic two-form in 



n = -i [dz{d)dz + dz(a)dC + d(,{^)dz + dC{y)d^ , 



(30) 
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where 



It then follows immediately that 



(1 + |Z|2)2 '"^ l+lzP'"^ l + |z 



d/i ^ dfj, (X sdet\\D,^ij\\dzdzdS,d^ cc — . ^ dzdzd^dS,, 

[1 + \z\ ) 



which agrees with Eq. ( [28l) . Poisson brackets with respect to ()30p then give us 

{^,^}PB = l-m^l^ {^JcovjpB ^ 0, 

which clearly implies that the amplitudes ^, ^ do not represent spinless fermions independent of lattice spins: vacancies 
and lattice spins are correlated due to the NDO constraint. For example, the destruction of a vacancy necessarily 
results in the creation of a lattice spin. As we see below, a calculation of a purely fermionic correlater involves both 
the ^ and z variables in a nontrivial manner. 

If we ignored the factor (1 — in the measure we would end up with a partition function of a spinless fermion, Z = 
1 + e^*^''. The conventional fermionic amplitudes /, / obey the conventional rules, {f,f}pB = 1, {f, Scov}pb — 0, 
with the symplectic structure taking on the standard form, 

dzdz - 
^ ^ I I I2^2 ~ ■ 

(1 + |z|2)2 

The composite field (z,^) parameterizes a point on a supersphere, 5*212 ^ CP^^^. For the conventional spinless 
fermions coupled to su(2) spins the underlying phase space is instead given by a direct product of the ordinary 
sphere, 52 ~ CP^, and a complex Grassmann plane. 



IV. su{2\l) PATH INTEGRAL: PARTITION FUNCTION 

We turn now to a derivation of the path-integral representation of the partition function. A key ingredient in 
constructing the CS path integral is the resolution of unity in HpUys, 



dM|z,0(z,ei, (31) 

where the measure is given by (I^Hl) . This equation can be used repeatedly to compute (z, ^|e^^^|z, — ^), considering 
to be a local on-site Hamiltonian expressible in terms of the sm(2|1) generators. Following standard procedure, let 
us break up the interval [0, /3] into N small pieces of length e = P/N, N ^ 00. Then Eg. dT^ can be rewritten in the 
form 

r r ^ 

tre-^" ^ J dfij Y[dfik{z,m){N\N -l)---{0\z,-Oe-'^>^"^'''''-^^ +0{e^). (32) 



fe=0 



Here 



ff(fc,fc-l) = iffl^^, \k):=\zk,a), zk = z{ek),^k=a^k), k = l,...,N, (33) 
and equation ([TO)) tells us that 

J d^l{z,C\N){0\z,-O=tr\N)(0\ = {zo,^o\zn,~^n). 
Finally, integrating over d^ in ([5^ yields 

tre-^tf^ d^,ol[dfik{k\k-l){zo,UzN,-^N)e^'^^''^''^''-'\ (34) 

7, 1 



fc=i 
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From now on we drop the 0{e'^) contribution to the partition function having in mind that the continuum hmit will 
be taken eventually. We further notice that the kernel {z,£,\z' acts as a delta function with respect to measure d/i. 
To see this consider any vector that belongs to Hphys- Then resolution of unity implies that 



which in components reads simply 



Having this in mind, Eq.(j34|) becomes 



tre 



Uo— ^^iV: ^0 — — CiN 



k=l k=l 

Performing here a formal time-continuum limit e — >■ yields 



. N N 

I ]^d/ifcexp^[log(/c|fc-l)-ei/(A:,fc-l)] |^„=^„,5„=_j„ 



tre 



-I3H 



Ldr 

Jo 



where 



and 



L = -{z,i\— + H\z,0 



The SU {2\1) symplectic potential explicitly reads 



{z,^\-—\z,0 = t; 



9t' 



(35) 



(36) 



(37) 



(38) 



(39) 



and {z,^\H\z,^) = Hcov 

We now specify Hamiltonian to be that of the U = oo Hubbard model given by Eq.([2]) and make the change of 
variables ([ST]) . We then get 



N 



n-p-p dz,{ek)dz,{ek)di,{tk)di,{tk) ^ f , -.c r , \\ 



i fc=l 
N 



k=l 



2Tri{l + Zi{ek)z,{ek))^ 
J2^ogA,{k) ~ eH^ov{k,k ^ 1) 



(40) 



accompanied with the boundary conditions zg = zjy, = ~^n- Here 

H,ov (fc, fc - 1) = t ^ (efc)C. (e(fc " 1)) i^k)\z, {e{k - 1)) - ^^ i^k)(^ i<k~l)), (41) 



and 

Mk) = 



^ n .^M2^ (^ " \u^k)U^k)){l \u<k l))6(e(fc - 1))) + C.(efc)C.(e(fc - 1)). 
A/(l-|-|z,(efc)|2)(l + |z,(e(fc-l))|2) 2 2. 
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Let us recall that the indices i,j denote the lattice sites, whereas the index k numerates time slices only. 
In the continuum e — > limit, Ea. (HD|) becomes 

tre-P" = J Dfi{z, efo ^'^'^''^^ (42) 

where 

i>>.fee)=n '"-'1Sitt'p!''''^' "-'^^-' 

i.r 

stands for the measure with the boundary conditions, Zi{0) — Zi{p), ^i(O) = —^i{p). The Lagrangian now reads 

L ^ (^) - E {^r + ^^ + «af ^) ^^ - Hcov (44) 

i i 

The first piece of the action involves the time component of the Berry connection to be discussed later, 
In the time-discretized representation, it reads 



»QW(efc) = tak = ' , '^"f = \og{zk\zk-i} + 0{Sl), 5k = Zk-Zk-i. (45) 
2 l + \zk\^ 

The dynamical part of the action is given by 



Hcov = tY.^M^^l''^) - ^IIC^C^- (46) 

i^j i 



Here {zi\zj) stands for the product of the spin coherent states. 



1 + ZiZj 



(47) 



The covariant symbol of the on-site electron spin operator reduces to 

Qr = ^r(i-66). (48) 

Now we move on to a practical calculation of path integral (H^ . However, to be more accurate we do that right on 
the time lattice. As we see below, that will enables us to clarify some subtle points concerning the structure of the 
path- integral action which is not seen in the naive continuum limit. 

The major problem is of course the factor Hfcll ~ CfcCfc) presented in the on-site measure. Since one can rewrite 
this factor as 



exp ^- ^ Cfcffc j ^ cxp X! ^kS.ki^ exp ^^dr 



/3/N -> 0, iV -> oo 



it might seem that in the continuum limit it simply amounts to an additive renormalization of the chemical potential 
and as such can be discarded |14||. However, this is not the case. To start with, the renormalization is infinite 
which makes the whole procedure rather formal. Besides, in the chemical potential term from representation (1411) the 
arguments of the f and ^ fields are shifted by one time step, whereas the factor in the measure involves the ^ fields 
at coinciding time moments. This slight difference happens to affect the final results drastically. 
To see this, consider path integral (|42|) with the action 

/■/3 rl3 _ , . 

S = J ia^^\T)dT- j C(9r + /i + ia^"'jCrf-r- (49) 
With this action the path-integral evaluation of the partition function should reproduce the earlier result 
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which we aheady calculated in Eq. ([29)) using ordinary integrals. By a L'^(l) phase transformation of the fermionic 
fields the potential a^^') (r) can be brought into a time independent form, 



where 



Note that 0(0) = (/'(/3). The effective action then becomes 



(50) 



S 



(51) 



where ^ = + ia^^^dr. Here we cannot simply integrate out the ^ fields in (I5ip by standard means, since now 

the measure contains the extra factor, exp (— X]fe'?fcCfe)- To figure out how the path integral works in this case we 
must look back at the defining time-discretized representation PO]) . It becomes 



tre 



N 

n 2^z(l + ZfeZfc)2 



dzkdzkd^kd^ke""' _c 
■e ^, 



where 



N 



k=2 



(Cfc - Cfc-l) , 1 , c 

^ l-Kk-1 +e4i 



fJ'i.N 



N 



(52) 



(53) 



fc=i 



Here we have taken into account that — ~^n- The last term in the action represents the contribution that originates 
from the extra measure factor, so that the integral over the fermionic fields reduces simply to 



N 

lim / f] dfkd£ke'^"'=^^'''^'"^' ^ lim det 5. 

•' k=l 



In our case the N x N matrix S reads 



S = 



/ 2 

-b 2 







-6 2 



where 6 = 1 — -^^. The determinant is evaluated to yield 



lim det S = lim 

N->oo N->oo 



(1 - ^m)^ 



(54) 



(55) 



If, in contrast, we absorb the extra factor produced by the measure in the chemical potential term, we, instead, end 
up with 



lim 

N—foo 



1 + 



- 1 + e 



which yields the familiar result for non-interacting spinless fermions with the chemical potential 



yU — > ^ = ^ + 
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The 2^ factor in Eg. ([55)1 is produced by the modification of the measure. As discussed earlier in our analysis of the 
time discretized representation, it cannot be absorbed in the chemical potential term. Besides, integrating out the 
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Fermi fields formally results in a divergent expression in the limit N — oo. However, as we show below the remaining 
integral in (j52p over the bosonic fields z^, precisely cancels such a divergence. 
Explicitly, we get 

TV _ ■ A" _ ■ 

tre-^^^"°= hm / FT det 5 = lim / TT ^^^^^^(2^ + e-^^e^i:. -) 

N^oo J A^i^ 27ri(l + ZkZk)^ N^oo J 27ri(l + ZkZkf 



N 



^Ll 27ri(l + ZkZkf 



Um 2^ /fr , -t - .. e-+e 



smce 



27ri(l + zz) 



= 1. (56) 



Using the result 

f 1^ dzkdzk _ fyr dzkdzk 

J JJ- 27rz(l + ZkZk)^ J \^ 27ri(l 



fc=i 



27rz(l + ZkZk)^ 



^k\Zk-l) \zo=zm 



^ dzkdz, 



27rz(l + ZkZkY 



{2s + l)dzdz f 2dzdz 



where 

/ 27ri(l + zz)2 '""■■^"i"' " 7 27ni(l + zz) 
finally gives 

tre-^^^™ =2 + e^^^ 

as desired. On the other hand, if we suppress the extra measure factor we find instead the incorrect result. 

At this point it is instructive to compare the action (|49p with its formal analog for conventional spinless fermions 
coupled to SU{2) spins through the Berry's potential: 

S,/f = J ia^°\T)dT-J f {dr + lJi + ia^°'^j fdr. (59) 
The standard path-integral representation of the pertinent partition function reads 

Zs/f - / i?A^f~e^=/^ (60) 

where 

M~ - n L"^i"+ff "^"(")^/(") (61) 

is a standard measure in the CS space generated by the basis 1^)5^1/2 I/)- Here |/) stands for the normalized CS 
associated with the fermionic algebra, 

V -L + 7/ V + // 
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where {f,P} = 1- Within Shankar's approach, /''' represents a hole creation operator. Accordingly, the one-particle 
state 11)^- corresponds to a hole excitation, whereas \0)p stands for a hole- vacuum state. However, this vacuum state 
is unphysical, since it is not present in the 3D on-site reduced Hilbert space for strongly correlated electrons which is 
spanned by the spin-up, spin-down and hole (vacancy) state vectors. 

The integral (l60l) can be calculated on the time lattice along the lines depicted earlier on and yields the divergent 
result. 



s/f 



N 



lim 2(1 + 2 



(63) 



This clearly indicates that action (|59| cannot represent a physical system. The Berry's potential in the last term of 
Eq. ([59l) takes into account the fact that the su{2) spins and spinless fermions are correlated within Shankar's theory 
due to the NDO constraint, even in the absence of a direct interaction between them. However, that correlation is 
not taken into a full account there. The projection of the spin-fermion Hamiltonian onto the physically constrained 
Hilbert space manifests itself both in the appearance of the Berry's phase in the fermionic action as well as in the 
modification of the measure in the path integral ()42|) for the partition function. 

It is also important to realize that the spin part of the measure in Eg. (1611) contains the extra factor of 2 = (2s+l)g^j/2 
in the numerator compared to the spin part of the on-site version of Eq. (j43|) . It comes from the normalized SU{2) 
spin measure as given by the resolution of the identity in the spin space (I58[) . The fermionic part of the measure (j61|) 
is standard and comes from the resolution of the identity in the fermionic Hilbert space, 



ermioji 



dfdf\f){f\ 



This factorization of the path-integral measure (j6ip into pure spin and fermionic parts reflects the fact that the 
Hilbert space of the whole system is represented by a direct product of the spin and fermion subspaces. In contrast, 
the measure (|43p comes from the resolution of the unity in the whole superspace CP^'^ as given by Ea. (|3T|) . It cannot 
be factorized into the spin-fermion parts and the resulting effective action represents a unique composed object. For 
example, the integration over just the fermionic amplitudes in ()52|) diverges. 
Finally, let us evaluate the Green's function 



(64) 



where 



-m 



(65) 



and 



Sn = 



(0) 



(66) 



Making the change of variables, ^(r) ^{T)e "^(■^), where 0(t) is given by ([SU)) . this can be brought into the form 



-i:0(T,)+i0(T,)g5o(2,C) 



(67) 



where Tq = j^ct, ce is integer, and Sq is given by (I5ip . This expression is a continuum limit of the time-discretized 
representation 



/• ^g, piafc-i0(T,) + i</)(rr) 



(68) 



where the inverse of S is 



S^^^ = (det5)-i 



det5 = 2^ + 6^. 



(69) 
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Hence, for q > r 



and this gives 



C-l ^ f !_ (70) 

9^ 2^ + 6^ ' 



AT^ooJ 27ri(l +ZfcZfc)2 

n(^.kfe-i>(n(^^i^fc-i))"'- (71) 

k—l k—r 



X 



Since 



N , N q 

fc = l ^ I K 



we finally get 



- ^r) e-'^(^''-^'-)(l - 5), (72) 

where the vacancy occupation probability 5 is determined by 

5 = (XOO) ^ 



2e/5'^ + 1 ' 
Similarly, for q < r, we find 

G(")(T,-r,) ~e-^(^'"^-)5. (73) 

Combining these results, we finally get 

G(°)(r - t') = e-^(^-^') [6I(t - t' - 77)(1 - S) - 0{t' - t + r;)(5] , (74) 

where the infinitesimal parameter 77 = 0'*' is a reminder that the variables ^(r') and ^(r) in the path integral (|42l) are 
associated with variables displaced by one time step, ^{ek) and 5(e(fc — 1)), respectively. At equal times we find 

-G(°)(0-) = (e(r)e(r)) = (X°")=5, (75) 

as it should. In contrast, for the conventional spinless fermions governed by the Hamiltonian H = jif'^ f one obtains 
Eg. ((74)) with 5 replaced hy np- 

Because of the extra factor nr(l ~ C(''')C (''')) the measure (j43|) . the fields ^, ^ separately do not have a direct 
physical meaning. They don't represent spinless fermions, for instance. However, the bilinear combination S,{t)£^{t') 
represents, at equal times, a physical observable, the covariant symbol of the vacancy number operator X^^ . We see 
that the spin degrees of freedom are nontrivially involved in the evaluation of a purely fermionic vacancy propagator 
([74)) . Notice, however, that this spin-charge correlation due to the NDO constraint does not prevent the spin-charge 
separation observed in the \D Hubbard model. The physical elementary excitations that represent separately the 
spin and charge degrees of freedom in the ID Hubbard model are not simply the z and ^ field excitations, but are, 
instead, nonlocal string objects composed, simultaneously, of both types of bare elementary excitations [T5| . 

Finally, let us say a few words concerning the meaning of the continuum path-integral representation (j42p . Strictly 
speaking, Eq.(|42|) is only symbolic: the measure in that path integral cannot be defined in a mathematically rigorous 
way. This observation has a deep physical meaning. If a continuum path integral could be defined as a bona-fide 
integral with respect to a genuine measure insensitive to discrete approximations, it would immediately provide a one- 
to-one correspondence between classical and quantum physics thereby making the latter an unnecessary ingredient. 
However, classical dynamics is known to give rise to different quantum theories if one follows different quantization 
schemes. A specific quantization scheme is encoded into a discrete representation of the continuum path integral that 
should be taken as its true definition. In our case the basic discrete representation is given by our Eq. (j40p . Since 

^{r + r^)^ar) + '^V + 0{rf ), ry ^ 0, (76) 
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it might seem that one could safely ignore that shift in the low-energy limit, ^ cx w — >■ 0, where w stands for a 
characteristic frequency of the fermionic degrees of freedom [l^. However, this is not the case: the path- integral 
variables are not in general smooth functions. There may exist different discrete approximations to one and the same 
continuum action that result in completely different low-energy dynamics. Some explicit examples can be found in 
p7| . Even in cases where the continuum limit of a path integral does make sense, i.e. in a semiclassical or perturbation 
theory, the relative shift of the arguments of the path-integral variables still cannot be ignored. The discontinuity of 
the correlators such as < £,{t)(,{t') >, at equal time arguments, should be dealt with according to rules following from 
the defining discrete approximation. As soon as the discrete approximation is fixed, no "operator ordering problem" 
shows up. In our case the order is fixed from the very beginning and manifests itself in the Hamiltonian function 
defined by Eq.iP]). 



V. ID EXAMPLE 



The U = oo Hubbard model ^ is known to be exactly solvable in \D. We show below that the exact ground-state 
energy can be recovered within the path- integral representation (|42j|44p . To start with, the ground state of the ID 
U = oo Hubbard model is known to be degenerate with respect to spin. To calculate the path integral we can therefore 
choose any spin configuration, e.g., the ferromagnetic (FM) one. We thus put Zi = Zj in Eq. (P^ . which yields 

Zu=oo^ I D^i{z,S,) e^^^^^\ (77) 



where the measure is given by Ea. (|^5)) and the action reads 

S = J tria'^^Ur-J x [dr + ia'^"^ + Tj xdr. (78) 

Here T and a^^"^ are the Ng x Ng matrices, 

{a^°%=afS.„ (T),, =t,,-M%, (79) 

with i,j = l,2,..^Ns numbering the lattice sites. The vector x = i^iy£.2, ^n^Y ■, and the trace is taken over the 
lattice site indices. On the time lattice, the integral over the x fields can be carried out to yield 

Z._= hm / n n f v'^f-?.wC det(2^ + e-^-~^^.^n 



Since the matrix ||2^e*Sfc°'° -t- e^^|| at = is diagonal, Ea. (l5Ul) reduces to 

Zu=oo{t = 0) = (2 + e^^)^^ = e^^^=(l + o(l)), /3 ^ oo, ^ > 0, 

where Ns is fixed. This result is a direct consequence of Eq. ([57)) . Accordingly, at ty ^ the negative eigenvalues of 
the matrix T determine the asymptotic behaviour of the partition function. Consequently, 

Zu=oo^ n e-'^^''(l+o(l))> /3^oo, 

p:Tp<0 

where Tp — —tp — /i, tp = 2tcosp, p G BZ. The ground-state energy then becomes 

Egr/N, = ~—sin{TTS), ^ = -arccos(-^), i > 0, (81) 
which coincides with the exact ID result for the U = oo Hubbard model UM. 
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VI. DOPED ANTIFERROMAGNET 



Let us now turn our attention to a derivation of the low-energy effective action of a doped AF, starting riglit from 
the microscopic t — J model, 

Ht_j = Ht + Hj = -tY^ 45,, + J ^ • Q„ (82) 

ija ij 

where Qi stands for the local electron spin operators given by Eq.Q and J > describes the nn exchange interaction. 
The parameter J ^ 0{1/U) and the bare constants are chosen such that t ^ J. 

Because of the NDO constraint there are no charge fluctuations at half filling {S — 0), and precisely in this limit, 
the t — J model reduces to a Heisenberg AF model 

Ht-j ^ Hj=° ^jY,S^-Sj, (83) 

with no restriction on J, the sole energy scale, apart from its positive sign. In the low-energy long-wavelength limit 
Hj gives rise to the action of the nonlinear sigma-model. 

Proceeding as discussed in the preceding section we arrive at the representation of the t — J partition function, 

Zt-j = I Dfi{z,0 er^*-^(^'«)''^ (84) 

where the measure factor Dfj,{z^^) is given by Eq. (j43l) . The Lagrangian now reads 

Lt.j = Yl ^^i^V) - E {^r + /i + ^af^) 6 - (85) 

i i 

where 

The summation in (1861) is extended over nn and nnn sites. As a result wc find 



= (z(rl)|z(rl + 5f)) = 1 + i^z{n)\^\z{n)W^O{8T») = 1 + (z(rO||| A + ^^\,(^r^^))Sf + 0{Sr^) 

= 1 - iaidf+ 0{Sr^)^ exp(-Wi(5f) +0{d7^). 

where fj = fi + Sr. Vector a represents the spatial components of the pull-back of the Berry's connection one-form, 
ip*A = a, with if : (r, f) — )■ z(r, f) and 

A = i{z\d\z). 

Here d stands for the exterior derivative, 

a = dz— + dz — . 
oz oz 

The effective hole action becomes 

St ^Y r*«r^W'^^-E (a.+/i + mf^)e.-dr 

' ^ f H '^f.+Sr<f.e-^'''''dr + 0{57»). (87) 

JO - X- 

Here (r) := a^{fi,T), v = 0,a;, j/, with ia^ = —{zi\-^\zi). Note also that Zi := z{ri) and rj = + where Sfcca 
with a being the lattice spacing. 
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The whole action should be complemented by the J-term. Right at half filling the only surviving term gives rise to 
the nonlinear sigma-model action to describe the quantum AF. The perfect Neel ordering prevents interlattice hole 
hopping. This follows directly from Eq. (1551) . Namely, the AF long-range order, Si — —Sj — S, implies 

Zi ^ z, Zj ^ -1/z, ieAjeB. (88) 

In view of (|47|. this results in {zi(zA\zj^B) ~ {z\^l/z) — 0, and there is no hopping between the A and B sublattices. 
Only intralattice nnn hopping of vacancies is possible in this case. One can also check that 

{z\d\z) = -{-l/z\d\-l/z), (89) 

which means that aj'g^ = '^o.ieB- 

Suppose now that we lightly dope the AF with holes, with the dominant energy scale in the problem continuing to 
be the exchange coupling J. We assume that the AF ordered lattice spins are slightly perturbed by a small amount 
of vacancies. In this case for a small enough hole concentration S, one gets 

z,^-l/zj+0{d), 

where i,j are the nn sites. This in turn implies 

{z,\z,) ^0(6), 6^0. 

Therefore , very close to half filling, the interlattice hopping effectively results in the renormalization of the hopping 
amplitude t ^ St and, as a result, this term can be discarded. 

Summing this all up, the partition function of the t — J model close to half filling reads 



Here 



Tie A Tie A 



J J2 s^rsf:isf, j = j{i- (91) 



and the measure is given by Eq. 



Within Shankar's approach, the underdoped AF is described by the partition function 



Shankar ^ j g/o (92) 



Here 



TiEA Ti^A 

- t J2 k+Sfhe-'^-'"' +iA-^B,a,-^- 



J ^ ^ Sft^ Sp._^.gp, (93) 



ri ,(5r 

and the unconstrained spin-fermion measure reduces to 



7-, fermion TT 2 d^:^ (t)^^^ (t) - , , ,j, , ^ /nA\ 



i(l + ZiZi) 
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Here the amplitudes denoted by fi describe the conventional spinless fermions, whereas the Zi fields correspond to 
the su(2) spins. This is Shankar's result for the underdoped AF (with s = 1/2). 

Formally, the Lagragians (|9T|) and (|93| appear to be identical. However, due to the difi^erent measures, the cor- 
responding partition hmctions are completely different from each other. The fermionic extra term in our measure 
(|43| cannot be ignored. Physically, it reflects a rearrangement of the underlying Hilbert space induced by the NDO 
constraint. This clearly cannot be treated perturbatively. In a true spin-fermion Hamiltonian, the on-site Hilbert 
space is represented by a AD direct product of the 2D spin and the 2D spinless fermion subspaces. However, the 
NDO constraint reduces it to a 3D on-site Hilbert space composed of three state vectors: the spin-up, spin-down and 
the vacancy states. As shown above, it is precisely the extra measure term that takes explicit care of that distinction. 

As an example, consider partition functions ([TO]) and (IMI) in the limiting case t — J = 0: 



Lt-At = J = 0) = J2 »a! (^) - E [dr + /i + *af ^) , (95) 

and 

^Shankar^^ = J = 0) = ^af\T) - ^ + ^Ji + laf^) M- (96) 

In view of Eqs.(| 



Zt^j{t = J = 0) = {2 + e-^^)''% 

as it should. On the other hand, Ea.([55| tells us that ^•S/ianfcar^^ = J = 0) diverges, which is inappropriate for a 
system with a finite number of degrees of freedom. If one, however, drops the af''^ potential in the fermionic part of 
Ea.(|96p. one arrives at the conventional partition function 

2^-'(l + e-'^^)^= 

which describes entirely independent uncorrelated spin and fermionic degrees of freedom. 

Only right at half filling, 6 — 0, our representation (^0]) and Shankar's theory become identical. This can be seen as 
follows. Since there are no holes in this limit, the fermionic amplitudes (/, /) describing their propagation throughout 
the lattice vanish identically. The partition function then becomes 

Sharikar ^ f L^^''-'='"-(5=0)rfr 



^5=0 = / i^Msn(2)eJo , (97) 

where 

j^Shankar^^ = 0) = ^ ^ (r) - iJ™" = 5™"^™", (98) 

i i^j 

and the normalized spin measure is given by 

i^M.„(2)=ll . (99) 

On the other hand, at i = 0, the projected t — J Hamiltonian reads 

ij i 

The half filling limit in this representation can be enforced by sending the chemical potential /j, to +00. This has an 
immediate effect in allowing only the zero eigenvalues of the local vacancy number operator X'j^'^ to survive. In this 
way the partition function becomes 

Zj^ J Df,{z,OeIo^'''\ (101) 

with 



jYsrS''°\ (102) 
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The measure DfjL{z,^) is given by Ea. (l43l) . Note that it is normahzed by the condition (|3T|) . so that the factor of 
2 = (2s + l)s=i/2 is missing in the numerator of its spin part. 

Notice that we cannot simply put ^ = ^ = in Ea. (|102p . This is because those fields have no physical meaning and 
do not represent holes directly as opposed to the amplitudes /, /. We must instead integrate them out in the limit 
^ — T' +00. We do that on the time lattice, which gives 

Z, ^ hm 2--= / n ft f-l^^t-^w'i^^ E - 1)]. (103) 

N-^co J 27rz(l + Zfe(j)zfc(j))^ ^ 

Here 



Hcov{k, fc - 1) = 81°'' [zk, Zk-i)Sj°''{zk,Zk-i). 

Note that the factor of 2^^= comes precisely from the fermionic extra factor in the measure (|43)) . This factor can 
be absorbed back into the spin measure to turn it into a conventional normalized spin measure (I99p . As a result, 
the integral over the spin fields z,z in (I103P becomes identical to that given by Eq. ([TT)) . The distinction between 
the representations (l90l) and (|92|) disappears at half filling since the corresponding on-site spin Hilbert spaces become 
identical. 

Finally, the Berry's gauge potential a"^ = {(p*A)'^ appears in the path-integral action (PT|) as an external gauge field 
and has no dynamical role. The gauge-theory approaches to treat the t — J model basically fall into two categories. 
The first one emerges from the slave-particle representations of the constrained electron operators. The idea of 
that approach is to re-express these constrained operators in terms of the standard boson/fermion bilinears. This 
is equivalent to the so-called oscillatory representations of su(2|l) superalgebra. For example, the physical electron 
operator can be represented by a product of a spinful boson and slave spinless fermion 

with standard commutation/anticommutations rules. However, this representation clearly increases the on-site number 
of degrees of freedom by a factor of 2. The emergent U{1) local gauge field, fi, bia — ?> e*^*/i, e'^'&io-, takes care of one 
redundant degrees of freedom (by fixing a gauge), while the NDO constraint 

a 

takes care of the other. However, the elementary excitations of the slave-particle fields do not represent physical 
excitations since they are gauge dependent. This compact dynamical U{1) gauge field is generated by the NDO 
constraint. Within our approach the constraint is resolved explicitly and there is no need for an imposing the U{1) 
gauge field theory. In fact, an explicit resolution of the NDO constraint within the slave-fermion theory results exactly 
in the su{2\l) path-integral representation 

The gauge potential = {ip* AiY which appears in Eq. ipTj) has a different nature altogether. It is driven by the 
fluctuations of the spin background present in any su{2) path-integral action regardless of the chosen coordinates and 
it takes care of the U{1) redundancy of the spin quantum state \z). From the geometric viewpoint, this can be stated 
as follows [19|]. The su{2) CS's can be considered as sections of the principle U{1) bundle P{CP^, C/(l)) frequently 
referred to as a magnetic monopole bundle. The base space of this bundle, M = CP^, appears as a classical phase 
space of the spin, whereas its covariantly constant sections \z) form a quantum Hilbert space for the spin, with 

V|z)=0, \/-~d + iA. 

Since the manifold CP^ ~ S"^ is topologically nontrivial, the monopole bundle is nontrivial as well. Notice that 
the one-form A does not exist globally, and any two locally defined gauge potentials are related by a U{1) gauge 
transformation. Fixing the spin CS by Eq. ([7]) amounts to fixing a local section of the bundle. 
Under the global canonical SU{2) transformations acting in the base space we get 

" 9=(\l)^Sm- (104) 



-VZi +u \ -V u 

When lifted to the bundle (quantum) space, this gives 



g^Ug: Ug\z.)^e-<^\z!), Q^-ilogJ (105) 
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The Berry's connection transforms according to 

A^A, + dQ. (106) 

and the Fermi field becomes 

^^ ^ e^'^M (107) 
leaving the whole action in Eq. ([90)) globally SU{2) invariant. 

VII. DISCUSSION 

In this section we discuss some physical consequences that follow from representations (|90|) and (|92l) . Let us start 
with Shankar's theory (|92p . This theory naturally allows for the hopping of the conventional fermions throughout 
the lattice in the absence of the local spin degrees of freedom. The corresponding partition function follows from the 
representation (j92[) if one discards the spin degrees of freedom. In the momentum representation, it reads 



ZSHar^kar ^ | ^ dUt)DU{t) exp j - f .U<dr ~ ep)fpdT j . (108) 



Here — t^; — /i, and tjj — 2t cos(p a), where a is a lattice vector and p G BZ. This path integral can be easily 
computed to yield a partition function for the conventional spinless fermions, 

^Shankar _ p^*^?) 
P 

Having this in mind, the low-energy long-wavelength limit can be taken explicitly to reduce Shankar's action (|93|) 
in ID to that of Dirac fermions coupled to the nonlinear sigma model via a compact U{1) gauge field [1]. Following 
the usual procedure for ID systems to take into account the low-energy fermionic excitations, we linearize the theory 
near the Fermi points dzkp. The spinless Fermi amplitudes ip = ] j \fa can be written as follows: 

^(n) = e^''^^'Vi(") + e-'^'''^^2{n). (109) 
Here index n stands for the lattice sites. The resulting action in the continuum limit reads 

Zyiankar = ^ D^Di,D^X,^(2^[z,z)e'-^''\ 

Sp — J [-ipAi—d - ia)^A + ipB{—d + ia)^l;B]dxdT, 

Se = -^J dxdT{cd^zd^z + c'^zz)+ieW, c = Ja, (110) 

where Sg is the sigma model action including the topological O-term and a — a'^^^- The chemical potential /i is 
incorporated in the theory through the relation kp — arccos(— /i/2t). The Euclidean 2x2 gamma matrices can be 
taken in the form 70 = ay, 71 = so that 75 = 17071 = <Tz- Notice also that ip = (tpi, V'2)*, "0 = V'Wo- 

An interesting observation concerning the representation ()110p is that the dependence on the parameter actually 
drops out from the theory [l[. This can be seen as follows. The action to describe the massless fermions coupled with 
the spin background fields z, z, 



Sp = j [ip(—d ± ia)'ip]dxdT, 

is invariant under a chiral U{1) transformation, 

iP e^^^'^ip, i) VSe^T^"^, (111) 
with (f> being a parameter. However, this is no longer the case for the partition function 

Zp = J D^pD^l^e^". (112) 
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The fields z{x,t), z{x,t) map a compactified 2D plane {x,t) homeomorphic to a two sphere S'^ onto a spin phase 
space which is also a two sphere, 5^ — >■ S*^ It is known that such maps can be classified by the integers W which 
define the number of times the the "space-time" sphere covers the "spin" sphere. Explicitly, 

[ da. 

The remarkable result in quantum field theory tells us that, at nonzero values of W , the classical chiral symmetry 
cannot be promoted to quantum level. In fact, 

Zf ^ e'^^Zp, (113) 

under the transformation (|111[) . This manifests the presence of a quantum anomaly. In view of the well known index 
theorem, the winding number W is given by the difference of positive and negative chirality zero modes of the Dirac 
equation, 

W^n+-n^. (114) 

This implies that the Dirac operator D = —d ± ia has zero eigenvalues at 7^ 0, which kills the fermionic path 
integral in (|110p . The theory survives only if W = 0. In this case 9 then multiplies zero and nothing can depend on 
it. 

The direct consequence of that observation is that the difference between integer and half integer spins goes away 
and the low-energy excitations in the spin sector of the (jllOp model become massive. This is in agreement with the 
mean-field theory of the ID t — J model obtained within the slave-fermion framework |20']. However, in the exact 
excitation spectrum of the t — J model, both spin and charge excitations are gapless in \D, at any hole concentration 
[Tsl [2l| . Specifically, in the limit J <^ t the ID Hubbard model reduces to a squeezed Heisenberg chain with an 
enlarged lattice constant a~a/{l — 5) and a renormalized superexchange coupling 

J-J^^-^) [}- 2n{l-S) 

For small momentum the energy varies linearly with momentum, e{k) oc Jk, fc — 0, leading to a linear-temperature 
contribution to the low-temperature specific heat. The predicted linear term in the specific heat has been experimen- 
tally observed in a ID organic molecular solid [l^l- This system can be described in terms of the ID Hubbard model 
with a transfer integral of 2.1 x 10~^ eV and an effective Coulomb interaction of 0.17 eV. The action (|110p predicts 
instead an exponential fall off of the specific heat and hence does not capture the low-energy physics of ID strongly 
correlated electrons. 

On the other hand, within the representation (l90l) Eq. (|108[) is replaced with 

^Ko = / n Di,yiT)D^fiT) exp j - ^ ^firMr) " E " ^p)^p(^)dT , (115) 

where the above discussed shift of the arguments of the path-integral variables is implicitly understood. The time- 
lattice computation of this integral reduces to Eg. (1551) with the parameter JI replaced by fi. As a result, the partition 
function (|115|) is found to diverge. This divergency is of course readily compensated when the spin dynamics is 
turned on. However, the fermionic degrees of freedom taken alone cannot be considered as well-defined physical 
entities independently from spin variables. They act as auxiliary degrees of freedom to describe a composite object - 
a constrained electron. This also means that we cannot use the standard theory to describe the low-lying fermionic 
excitation as given by Eas. (jl08p . (|109p . The arguments that lead to Eq. (|113l) are no longer applicable. Accordingly, 
the gapless low-lying spin excitations cannot be ruled out in the present case. 

We have not yet discovered how to perform the path integral (|90p in the continuum limit in an analytically 
trustworthily way. Integration over the ^ or 2: variables separately from each other does not in general result in a 
physically relevant effective action. For example, the action ()49|) corresponds to the constrained on-site Hamiltonian 
(im. Integrating out the fermionic degrees of freedom results in the divergent expression given by Eg. ([55)) . A strongly 
correlated electron system at finite doping can hardly be represented solely in terms of the su{2) spin operators. This 
agrees with the observation that the spin-spin correlators in the ID t — J model involve the conventional su{2) spin 



operators necessarily modified by fermionic "string" operators. Weng et al. [15[ have shown that the effects due to 
the squeezing and rearrangement of spin configurations in ID t — J model can be included by introducing a nonlocal 
fermionic "string" field. These processes of first squeezing the t — J chain and then rearranging the spin configuration 
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cannot be described perturbatively and must be taken into account by introducing string-like fields. Formally, this 
can be thought of as a manifestation of the fact that the constrained electron (Hubbard) operator cannot be split 
into a convolution of conventional fermion and su{2) spin operators. The bare charge and spin degrees of freedom 
are not independent and interact with each other very strongly due to the NDO constraint. Within the path-integral 
approach, this statement reasserts itself in the appearance of the extra factor in the measure of the path integral (|90l) . 

VIII. CONCLUSION 

Let us now summarize the main results of our work. If one assumes on some phenomenological grounds that 
the low-energy physics of a strongly correlated electron system allows for a description in terms of the conventional 
fermions coupled to the conventional lattice spins, one quite naturally ends up with the action first suggested by 
Shankar. Strong correlations are encoded there as follows. The theory comprises the spinless fermions which, even in 
the absence of the direct spin-fermion interaction, are coupled to the lattice spins through the Berry's phase potential 
ao as exhibited in Ea. (l93p . That interaction brings in a correlation between the fermionic and spin degrees of freedom 
due to the NDO constraint. However, this doesn't seem to be enough to account for strong correlations to full extent. 
The example given by Eg. dMl) tells us that this approach is not truly self-consistent. 

Our derivation of the appropriate low-energy action directly from the microscopical t — J model brings out an 
important deviation from the Shankar's theory. Although the corresponding actions look formally identical, the 
measures in the path integral representations for the corresponding partition functions as well as the nature of the 
germane fields are essentially different. The path-integral representation of the t — J partition function ((90)) involves 
the measure explicitly modified due to the NDO constraint. That modification ensures that we work in the physical 
reduced Hilbert space. As a result, the fermionic and bosonic amplitudes in the action (PT|) no longer correspond to 
the true fermion and lattice spin degrees of freedom, but rather represent a unique composed object that describes 
the constraint electron as a whole. Consequently, the doped t — J model does not admit a representation in terms of 
spinless fermions coupled to the local AF fluctuations via a compact U{1) gauge field even in the low-energy limit. 
The explicit resolution of the NDO constraint results in a rather involved path-integral representation for the t — J 
partition function (|M|) . It is not yet evident how one can proceed with its direct calculation, except in some trivial 
cases. This appears as an evident consequence of the strong coupling nature of the physical low-energy excitations of 
a system of the constrained electrons. 

Within the slave-particle approach, this problem is reflected in the lack of a controlled treatment for the emergent 
gauge field that strongly couples holons to spinons. The NDO constraint which actually gives rise to that gauge 
field [23j is of a crucial importance right in the underdoped ((5^1) region which, by no accident, is the most 
interesting region of the phase diagram. It therefore brings in intractable strong interactions in this region, which in 
turn makes analytical calculations completely uncontrolled. Essentially, the gauge theory based on the slave-particle 
representation examines fluctuations around a mean-field solution that corresponds to a mean-field treatment of the 
NDO constraint. In general, such a theory is not stable against quantum fluctuations and the self-energy corrections 
are infinite due to the low-frequency gauge field fluctuations. The slave particles are not truly present in the physical 
spectrum, and cannot be treated as quasiparticles weakly coupled to the gauge field Q . 

On the other hand, the Hartree-Fock approximations in the fermionic path integral for the Hubbard model, 

Htu = 51 ^-'i'y^'J" + ^ XI "^^^4' (116) 

are known to recover successfully some well-known mean- field approximations such as, e.g., Stoner magnetism {2^. 
Those approaches are essentially based on the Hubbard- Stratonovich decomposition of the [/-quartic term which 
preserves the global spin rotation symmetry [2^. The Hubbard-Stratonovich fields are then considered within the 
mean-field theory, which basically amounts to a saddle-point approximation accompanied by an integration over 
the corresponding Gaussian fiuctuations around its saddle point solution. This is equivalent to a random-phase 
approximation (RPA) and it is normally controlled by an appropriately chosen 1/iV/ expansion, where Nf stands 
for the number of the relevant field components. However, it is very unlikely that the hole dynamics in the strongly 
correlated regime, close to half filling, admits such a mean-field description. In fact, the physics of Nagaoka's phase 
can hardly be recovered within such a mean-field theory [13] . 

In our view, the most appropriate way to proceed is to address the problem of the low-energy dynamics of the t — J 
model directly in terms of the superfield {S,,z) S CP^'^ as dictated by Eqs. (|36ll39p . After all, it is this composite 
field that represents the true physical degrees of freedom - the constraint electron excitations. In this way, one is 
supposed to end up with a sigma model with the CP^'^ target space. In fact, the 2D nonlinear sigma models with 
the (7P"+'"-i|" target superspaces have been discussed to describe percolations, polymers as well as some other 



22 



problems in statistical mechanics j28| . This program can be explicitly carried out for a supersymmetric t — J model 
that exhibits a global SU {2\1) invariance at J = 2t. Given a gromid state of the lattice t — J model, the low-energy 
dynamics of the fluctuations around that vacuum state is supposed to be controlled by the nonlinear CP^'^ sigma 
model. However, in a physically relevant case, t 3> J, the vacuum state is still unknown (in fact, there exist different 
competing vacuum states depending on the doping regime they are associated with). Therefore, it is still not clear 
what kind of low-energy action actually emerges in these cases. 

To summarize, there is still no complete analytical effective theory to describe the low-energy dynamics for the 
underdoped Mott insulator in general. Such a theory should account for a simultaneous existence of a few competing 
nontrivial features: short-range AF order, superconductivity, uniform and modulated spin-liquid regime and a pseu- 
dogap phase [29j. The important ingredient behind that picture is the local NDO constraint that essentially affects 
the low-energy physics close to half filling, and no trustworthy mean-field treatment is yet available to tackle this 
problem. 
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